第 38 卷 第 3 其 工程 热 物理 学 报 VoL38, No3 
2017 年 3 月 JOURNAL OF ENGINEERING THERMOPHYSICS Mar., 2017 


旋 技 射流 雾 化 的 直接 数值 


模拟 研究 
如 长 孝 罗 地 柳 建 人 


(浙江 大 学 能 源 清 洁 利用 国家 重点 实验 室 ， 杭 州 310027) 


摘 要 雾 化 现象 涉及 到 气 液 界 面 复杂 的 拓扑 变化 ,使 得 雾 化 过 程 极其 复杂 . 目前 针对 一 次 雾 化 的 机 理 研 究 及 模型 尚 不 
够 准确 ， 由 于 数值 计算 方法 的 快速 发 展 ， 赵 来 越 多 的 研究 者 采用 数值 方法 来 研究 一 次 雾 化 本 文采 用 大 规模 的 直接 数值 
模拟 手段 和 Level Set 追踪 气 液 界面 的 方法 模拟 旋转 射流 的 一 次 雾 化 过 程 , 以 揭示 旋 流 雾 化 的 机 理 和 界面 不 稳定 波 的 贡献 . 


关键 词 Level Set 方法 ; 
中 图 分 类 号 : O359 


一 次 雾 化 ; 旋转 雾 化 


文献 标识 码 : A 文章 编号 : 0253-231X(2017)03-0553-04 


Direct Numerical Simulation of Swirling Liquid Atomization 
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Abstract Liquid atomization exists in a wide range of applications, such as printing, agriculture 
and combustion devices. However, liquid atomization involves topology change of gas-liquid interface 
that makes it complex. By now, it is still not fully understood of the atomization mechanism and the 
atomization modelling is not accurate. With the development of computational capacity and interface 
capturing method, more and more researchers perform the numerical simulation of atomization. The 
aim of the present work is to study the primary atomization by using the Level Set interface cap- 
turing method to reveal the swirling atomization mechanism and the contribution of surface instability. 
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比如 混合 、 打 印 技术 、 食品 工程 、 农 业 、 制药 过 程 和 
发 动机 内 部 的 雾 化 燃烧 过 程 等 . 雾 化 过 程 极其 复杂 ， 
涉及 江 流 、 界 面 运动 、 液 体 结构 的 破碎 和 合并 等 过 
程 。 雾 化 过 程 分 为 一 次 雾 化 和 二 次 雾 化 ,一 次 雾 化 
为 液 柱 破碎 为 液 桥 和 液 滴 的 过 程 ， 二 次 雾 化 为 液 桥 
和 液 滴 进一步 破碎 为 更 小 的 液体 结构 的 过 程 。 虽然 
存在 一 些 雾 化 模型 ， 但 是 没有 一 个 统一 有 效 的 模型 
适用 于 整个 雾 化 过 程 和 模型 的 精度 有 限 ， 这 促使 我 
们 针对 雾 化 过 程 进行 研究 . 

在 雾 化 过 程 中 ,表面 不 稳定 波 的 发 展 发 挥 着 重 
要 的 作用 , 其 中 存在 两 种 主要 的 界面 不 稳定 性 , 即 开 
尔 文 - 效 姆 赫兹 (KH) 不 稳定 性 和 瑞 利 - 泰勒 (RT) 
不 稳定 性 。KH 不 稳定 性 发 生 于 气 液 界 面 存在 速度 
差 的 条 件 下 , RT 不 稳定 性 是 液 相 渗 入 气相 的 横向 不 
稳定 性 . Kim et al. 的 直射 流 模拟 结果 与 Marmot- 
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不 稳定 性 产生 了 轴 向 波动 而 RT 不 稳定 性 产生 了 周 
向 波动 ， 导 致 液 桥 的 形成 。Herrmann 加 通过 模拟 
研究 横 射 流 发 现 ，KH 不 稳定 性 导致 迎风 侧 袋 状 破 
碎 ，RT 不 稳定 性 导致 近 喷 口 区 域 液 柱 附 近 液 桥 的 形 
成 . 然而，Xiao et al. 四 发 现 液 柱 上 游 表 面 波 的 发 
展 是 由 于 RT 不 稳定 性 导致 的 , 这 和 Herrmann 的 结 
论 相 悖 . 此外，Shinjo and Umemura 回 对 于 横 射 流 
的 研究 表明 Tollmien-Schlichting (TS) 不 稳定 性 是 导 
致 液 桥 形成 的 主导 机 理 . 但 是 针对 旋转 射流 的 研究 
很 少 , 对 KH 和 RT 不 稳定 性 的 贡献 尚 不 清楚 , 并 且 
旋转 射流 对 于 燃气 轮机 和 内 燃 机 的 空气 燃料 混合 和 
火焰 的 稳定 有 着 重要 的 作用 ， 本 文 针对 该 设置 进行 
研究 . 

本 文 在 发 展 的 守恒 型 LS 方法 基础 上 , 对 旋转 射 
流 雾 化 进行 了 大 规模 的 直接 数值 模拟 ， 对 一 次 雾 化 
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过 程 进行 了 详细 的 分 析 ， 并且 揭 示 了 界面 不 稳定 波 
1 数值 方法 
1.1 不 可 压缩 Navier-Stokes 方程 

气相 和 液 相 的 不 可 压缩 Navier-Stokes 方程 如 


下 : 
Ou 1 1 。 
VY = i Vo [vu Vu ])+9 (1) 
Op 
B+ VP=0 (2) 


其 中 ，ww 为 速度 ，p 为 密度 ，p 为 压力 ，g 为 重 
力 加 速度 ，/ 为 动力 和 慕 度 。 气相 和 液 相 的 物性 参数 
均 为 常数 ， 即 气相 中 p = pg，/ = kg; 液 相 中 
Pp = 0 /= 内。 界面 处 ， 速 度 是 连续 的 ， 而 物性 
存在 阶 路 , 即 [or = pi 一 pg, [hr = pW 一 Kg。 压力 阶 
跃 表 示 为 


pr=ort2Hrn Vu.n (3) 


其 中 , o 为 表面 张力 系数 ，« 为 界面 曲率 , n 为 界面 
单位 法 向 量 . 

Navier-Stokes 方程 的 离散 基于 均匀 交错 网 格 ， 
其 中 压力 值 和 Level Set 标量 储存 于 网 格 中心 ， 速 
度 值 储 存 于 网 格 面 中 心 。 对 于 Navier-Stokes 方程 的 
离散 ,空间 离散 格式 采用 二 阶 精度 的 有 限 差分 格式 ， 
时 间 推 进 采 用 二 阶 精度 的 半 隐 式 克 兰 克 - 尼克 尔 森 
格式 . 
1.2 Level Set 方法 

本 文中 Level Set 方法 用 于 捕获 气 液 界 面 , Level 
Set 允 数 $ (zx,t) 的 零 等 值 面 即 为 气 液 界面 。 通 稍 
Level Set 函数 $ 定义 为 到 界面 的 符号 距离 函数 , 即 : 


I$ (2,0)| = 12 — zrl (4) 


其 中 ,zr 为 界面 上 到 x 最 短 距离 的 点 . $ > 0 的 区 
域 为 液 相 , 5 < 0 的 区 域 为 气相 。 
界面 单位 法 向 量 和 曲率 的 计算 非常 简便 ， 


V9 
"7 oq ® 
上 一 一 V .也 (6) 


然而 , 传统 的 Level Set 方法 会 有 严重 的 质量 损 
失 问 题 . 这 里 我 们 采用 Olsson 和 Kreiss [gj 双 曲 正切 
函数 来 代替 符号 距离 函数 ， 


其 中 ,se 为 正切 厚度 . 界面 的 输 运 由 Level Set 输 运 
方程 确定 ， 
Vy. (up) =0 (8) 
值得 注意 的 是 , 也 不 能 始终 保持 双 曲 正切 的 
形式 , 这 会 使 界面 发 生 非 物理 的 变化 和 质量 损失 . 因 
此 , 需要 引入 重新 初始 化 方程 , 这 需要 求解 如 下 方程 
使 其 达到 稳定 状态 ， 
vWD) = vv) 9) 
其 中 , 7 为 伪 时 间 . 每 一 步 迭 代 过 程 中 将 上 式 算 至 稳 
为 了 进一步 提高 原 有 LS 方法 的 质量 守恒 性 , 我 
们 采用 了 根据 界面 曲率 来 进行 质量 修正 的 方法 , 具 
体 步 又 详 见 门 . 
1.3 Ghost Fluid Method (GFM) 方法 
Navier-Stokes 方程 的 离散 中 一 个 重要 的 问题 是 
对 于 压力 梯度 的 离散 。 对 于 式 (3) 中 界面 处 压力 的 
不 连续 性 采用 GFM 方法 加 来 处 理 . 假设 界面 工 的 
位 置 zr 位 于 网 格 点 zi 和 zi+l 之 间 , 压力 泊 松 方程 
中 的 压力 离散 可 以 表示 为 


Ee 
Orz \0 07 


1 
pr (Plit1 > Peg,i) = 


gi 


) 

— (Pg,i — Pg,i— 

Pg - [plr 10 
Mp2 p*Ax? a0) 


其 中 , p* = pg 十 Pp1(1 一 ,$= (xr -72i)/AY. 
2 数值 模拟 

射流 形式 为 液体 从 环形 面 以 10 m/s 的 轴 向 速 
度 和 10 m/s 的 切 向 速度 喷 入 ， 环 形 区 域 厚度 为 0.1 
mm， 因 此 ， 喷 嘴 直 径 D 为 0.4 mm。 计 算 区 域 为 
10Dx10Dx10D, 网 格 数 为 10243. 入 口 和 出 口 分 别 
为 速度 入 口 和 对 流出 口 边界 条 件 。 气体 和 液体 的 物 
性 参数 为 pg= 40 kg:m ,pi= 800 kgm 3，Jg= 
1.6x10-4kg:m-l.s-!, I= 1.6x10-3kg.m-1.s-1, 0o = 
0.036 N.m-!。 两 个 独立 的 无 量 纲 参数 分 别 为 Re = 
piLV/m = 2000 和 We = pLV?/o = 222. 旋 流 雾 化 
的 模拟 结果 如 图 1 所 示 ， 可 以 看 出 旋转 液 膜 发 展 为 
液 桥 并 破碎 为 液 滴 的 级 串 过 程 . 
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0° 
图 1 旋 流 雾 化 气 液 界面 


Fig. 1 Interface of swirling atomization 


3 结果 与 讨论 
3.1 界面 结构 

在 射流 初期 , 在 项 部 前 沿 后 形成 回流 区 , 与 Shinjo 
et al， 回 模拟 类 似 ,， 在 芯 菇 状 液 膜 边 缘 脱 落 形成 一 
部 分 离散 的 液 滴 . 之 后 液 膜 向 下 游 发 展 ， 光 滑 的 液 
膜 保持 一 段 距 离 ， 随 后 在 周 向 产生 不 稳定 波 ， 不 稳 
定 波 的 发 展 促使 液 腊 沿 轴 向 拉 伸 成 许多 液 桥 ， 这些 
液 桥 进而 断裂 形成 液 滴 ， 即 液 膜 的 破碎 呈现 级 联 的 
特性 . 因此 , 液 滴 的 大 小 和 液 桥 的 直径 呈 相 关 性 . 在 
下 游 区 域 ， 液 滴 的 凝 并 和 破碎 发 生 ， 使 雾 化 更 加 复 
条. 在 液 膜 的 发 展 过 程 中 , 由 于 回流 区 的 影响 , 雾 化 
角度 随 着 时 间 逐 渐 增 加 ， 最 终 稳定 在 一 定 值 . 值得 
注意 的 是 ， 旋 流 液 膜 形成 的 回流 区 使 一 部 分 液 滴 向 
上 游 移 动 ， 进 而 撞击 液 膜 ， 致使 液 膜 二 次 破碎 . 

与 单 相 气相 流 场 相 比 ， 液 体 的 引入 增加 了 流 场 
的 复杂 性 . 当 z/D = 1 ~ 2 时 ， 轴 向 速度 已 接近 
零 . z/D = 4 时 , 出 现 回流 区 , 回流 区 的 出 现 是 旋 流 
流动 的 一 个 重要 特征 . 弥散 的 液 滴 尾部 增加 了 气体 
的 扰动 , 使 流 场 中 出 现 许 多 的 小 涡 结构 . 
3.2 液 桥 形成 和 破碎 

随 着 液 膜 向 下 游 发 展 ， 液 膜 的 厚度 逐渐 减 小 ， 
不 稳定 波 逐 渐 增 强 。 当 不 稳定 波 达 到 液 膜 厚度 量 级 
时 ， 液 膜 的 破碎 产生 . 在 旋 流 流动 中 ， 主 要 的 不 
稳定 性 包括 Kelvin-Helmholtz (KH) 不 稳定 性 和 离 
心 不 稳 定性 外. KH 不 稳定 性 是 由 气 液 界面 的 剪 切 
造成 的 ， 而 离心 不 稳定 性 是 由 周 向 运动 造成 的 . 对 
于 气 液 两 相 流 而 言 ， 还 应 考虑 由 于 密度 差异 引起 的 
Rayleigh-Taylor (RT) 不 稳定 性 . 在 该 模拟 中 , 沿 轴 向 
的 Kelvin-Helmholtz (KH) 不 稳定 波 并 未 明显 出 现 ， 
沿 周 向 的 离心 不 稳定 波 和 RT 不 稳定 波 在 液 膜 端 部 
出 现 ,导致 液 桥 的 形成 . KH 不 稳定 波 未 形成 可 能 与 
旋 流 数 和 其 本 身 并 不 显著 有 关 . 

进一步 , 采用 速度 梯度 张 量 的 第 二 不 变量 & 来 
表征 涡 结构 。 从 图 2 可 以 看 出 ,高 涡 量 集中 于 液 膜 


破碎 为 液 桥 和 液 桥 破碎 为 液 滴 的 阶段 . 在 这 个 阶段 ， 
涡 结 构 决 定 了 液 桥 的 形成 方向 。 可 以 看 到 与 流向 牌 
直 的 涡 结构 (图 2(a)) 产生 与 流向 垂直 的 液 桥 ， 这 些 
液 桥 形成 了 与 流向 平行 的 液 桥 的 连接 . 与 流向 平行 
的 涡 结构 (图 2(b)) 产生 与 流向 平行 的 液 桥 . 


图 2 界面 和 涡 的 关系 ( 深 灰 色 代 表 涡 量 @, 浅 灰 色 代 表 气 液 界 
面 ) (a) 垂直 涡 ; (b) 流向 涡 
Fig. 2 Interaction between interface and vorticity (deep grey 
and shallow grey represent vorticity @ and interface, 


respectively) (a) vertical vortex; (b) stream-wise vortex 
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图 3 液 桥 的 破碎 


Fig. 3 Breakup of ligaments 


在 液 桥 形成 以 后 ， 液 桥 逐 渐 被 拉 伸 为 长 条 状 结 
构 ， 如 图 3 所 示 ， 这 些 液 桥 在 表面 张力 、 离 心力 、 
空气 动力 和 黏 性 力 的 共同 作用 下 ， 形 成 螺旋 状 的 前 
进 方式 ， 产 生 形 同 Capillary 破碎 模式 ， 形 成 多 个 
液 滴 。 在 下 游 区 域 ， 液 滴 的 碰撞 合并 、 液 滴 的 分 离 
以 及 液 滴 撞击 液 膜 等 现象 的 出 现 ， 增加 了 雾 化 的 复 
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杂 性 . 通过 求解 式 (12)、(15) 和 (16)， 从 而 得 到 波 数 


3.3 界面 不 稳定 性 分 析 

当 液 膜 发 展 到 下 游 时 , 随 着 液 膜 厚度 的 减 小 , 不 
稳定 波 逐 渐 增强 。 当 不 稳定 波 的 幅 值 达到 液 膜 厚度 
尺度 时 , 液 膜 开始 破碎 . Gallaire 和 Chomaz [9| 指出 
KH 不 稳定 性 和 离心 不 稳定 性 是 旋 流 流动 形式 中 两 
种 主要 的 不 稳定 性 . 此 外 , 在 该 气 液 两 相 流 中 , 还 应 
考虑 由 于 密度 差异 引起 的 RT 不 稳定 性 . 如 图 1 所 
示 , 该 模拟 条 件 下 , 在 射流 方向 并 未 出 现 KH 不 稳定 
波 , 这 归 因 于 表面 张力 换 制 波动 的 产生 。 定性 来 看 ， 
离心 不 稳定 性 和 RT 不 稳定 性 是 导致 液 膜 破 碎 为 液 
桥 和 液 滴 的 主要 因素 . 

定量 方面 ， 我 们 采用 线性 不 稳定 分 析 的 方法 来 
分 析 由 于 RT 不 稳定 性 导致 的 液 桥 形成 过 程 ,其 
中 增长 率 可 以 表示 为 : 

os pl— Pg oe 
Pit pg 

其 中 , w 为 增长 速率 ，ae 为 液 相 到 气相 的 界面 加 速 
度 ,为 波 数 . 最 大 增长 率 时 的 波 数 可 以 通过 公式 
(11) 微分 得 到 ， 最 大 增长 率 时 的 波 数 也 是 最 容易 出 
现 的 ， 可 以 求 得 主导 波 数 为 : 
ac(pl — pg) 

30 
在 液 膜 破碎 区 液 桥 的 数目 N 和 主导 波 数 之 问 的 关系 
可 以 表示 为 : 

N 一 TD 入 = kwax De/2 

其 中 和 为 最 大 增长 率 时 的 波长 . 

对 于 KH 不 稳定 性 ， 表面 张力 和 惯性 力 是 主要 
的 竞争 作用 力 上， 波 数 可 以 表示 为 : 


(ea 


—k3 
pl 


ok 


pit pg 


(11) 


kmax,RT = (12) 


(13) 


w= BEAU — 
pl 


其 中 AU 为 气 液 界面 的 速度 差 . 主导 波 数 通过 求解 
dw/dk = 0 获得 : 


(14) 


kmax,KH i (15) 


对 于 离心 不 稳定 性 ， 由 于 加 速度 朝向 外 侧 从 而 增 大 
了 不 稳定 性 (9. 可 以 得 到 主导 的 周 向 波 数 为 : 


AR 一 4 其 中 dg 一 


CD 


9 和 257 分别 为 离心 加 速度 和 环 量 . 可 以 发 现 式 (16) 
与 式 (12) 类 似 . 


分 别 为 21.32, 47.14 和 536.68. 应 当 指 出 ， 有 RT 不 
稳定 性 得 到 的 波 数 与 实际 模拟 得 到 的 波 数 (24) 很 接 
近 ， 从 而 验证 了 RT 不 稳定 性 是 该 旋 流 工 况 条件 下 
的 液 膜 破碎 为 液 桥 的 主导 不 稳定 因素 . 图 1 中 周 向 
波 的 出 现 主要 是 RT 不 稳定 性 作用 的 结果 ， 这 和 文 
献 [1], [2] 中 关于 直射 流 的 结论 一 致 . 


4 结 论 


本 文采 用 大 规模 的 直接 数值 模拟 手段 和 Level 
Set 追踪 气 液 界面 的 方法 模拟 旋转 射流 的 一 次 雾 化 
过 程 ,结果 显示 RT 不 稳定 性 是 主导 液 膜 破 碎 的 主 
要 因素 . 
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